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1 . INTRODUCTION 

One  of  the  most  interesting  developments  in  the  recent 
history  of  containment  research  has  involved  the  realization 
of  the  importance  of  the  residual  stress  field  as  a signifi- 
cant feature  of  the  containment  of  an  underground  nuclear 
explosion.  This  stress  field  is  the  result  of  the  plastic 
deformation  of  the  geologic  material  during  the  formation  of 
the  explosive  cavity  and  its  subsequent  relcucation  to  cin 
equilibrium  state. 

Rimer showed  the  dependence  of  the  stress  field  on 
several  relevant  material  parameters,  and  his  results  have  been 
confirmed  and  extended  by  App,^^^  Terhune,^^^  and  Allen. 


All  investigators  have  shown  that  material  strength  is 
an  important  determinant  of  residual  stress  magnitude  for  a 
given  initial  condition.  Relatively  little  information  is 
available  about  the  strength  of  geological  materials  of  inter- 
est, and  such  information  is  difficult  to  obtain.  Therefore, 
an  interesting  question  involves  the  extent  to  which  strengths 
are  determined  by  other  measurable  material  properties.  Spec- 
ifically, if  it  could  be  shown  that  the  strengths  could  be 
estimated  reasonably  well  as  a function  of,  let  us  say,  water 
content,  sound  speed,  and  grain  density,  such  a relation 
could  significantly  extend  our  knowledge  of  material  strengths 
and  our  ability  to  estimate  residual  stress  magnitude  at  the 
site  of  future  contemplated  tests. 


An  effort  has  been  made  to  seek  such  a relation  by 
what  might  best  be  called  a brute-force,  least-squares  fitting 
approach.  No  preconceived  relationships  were  assumed  in  the 
process . 
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The  numerical  work  was  done  by  a least-square  fitting 
progreim  written  at  the  UCLA  medical  school  entitled  "BMD02R, 
Stepwise  Regression."  This  code  generates  least-squares  fits 
of  a set  of  dependent  variables  as  a function  of  a practically 
unlimited  set  of  possible  independent  parameters.  The  code 
examines  the  correlation  between  the  dependent  and  independent 
variables  and  chooses  the  most  significant  parameter  to  be 
used  in  generating  a least-squares  fit.  After  that  fit  is 
performed,  the  remaining  independent  parameters  are  again 
examined  and  the  second  most  significant  is  chosen,  and  a two- 
parameter  fit  is  generated.  This  sequence  is  continued  until 
further  fitting  becomes  statistically  insignificant.  During 
the  course  of  a fitting  cycle,  parameters  can  be  either  intro- 
duced or  removed  from  the  fit. 
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2.  TECHNICAL  DISCUSSION 


Material  properties  information,  determined  by 

TerraTek  for  area  12  tuff  samples,  was  used  exclusively  in 

this  investigation.  The  strength  parcuneter  used  is  the  stress 

difference  at  4 Kb  from  uniaxial  strain  tests  as  scaled  from 

r 5—91 

data  published  in  a series  of  recent  TerraTek  reports. 

These  strength  parameters  were  determined  for  the  48  drill 
holes  in  E,  N and  T tunnels  of  area  12,  listed  in  Table  I. 

A total  of  471  data  points  were  found  in  which  stress  difference 
at  4 Kb  as  well  as  a complete  suite  of  other  material  pro- 
perties information  was  available.  The  other  information 
included  as-received  and  grain  densities,  water  contents, 
and  longitudinal  and  shear-sound  speeds  determined  from 
ultrasonic  tests.  Porosity  and  saturation  values  were 
derived  from  these.  For  purposes  of  this  exercise,  porosity 
and  saturation  were  treated  as  independent  parameters  even 
though  they  were  derived  from  other  measurements. 

Figures  1-7  show  "strength"  plotted  as  a function  of  the 
seven  "measured"  quantities.  It  is  interesting  to  note  in 
every  case  that  the  relations  can  best  be  described  as  shot- 
gun patterns.  In  no  case  is  a clear  relationship  evident. 

Figure  1 shows  the  relation  between  strength  and  water  content. 
It  is  quite  true  that  the  field  of  data  can  be  limited  on  the 
upper  right  hand  side  by  a line  which  indicates  that  strength 
goes  to  zero  at  roughly  32  wt.  percent  of  water  in  tuff.  Such  a 

relationship  has  been  used  in  some  calculations  of  event 

parcuneters.  However,  that  is  only  one  limit  to  the  data  and 

can  heirdly  be  called  a representation  of  the  information 
available. 

This  set  of  471  strength  determinations  was  fit  by 
linear  least-squares  techniques  to  an  optimum  choice  from  113 
possible  parameters.  Specifically,  the  fit  used  the  measured 
parameters  themselves,  various  modulii  derived  from  them. 
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Sources  for  471  Data  Points  Used  in  a Study  of  the  Influence 
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OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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MASS  FRACTION  OF  WATER 


Figure  1. 
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OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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Figure  2. 
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OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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Figure  3. 
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OBSERVED  ROCK  STRENGTH 


OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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Figure  5 . 
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OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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Figure  6. 
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OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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Figure  7. 
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miscellaneous  ratios  of  these  modulii,  and  the  squares  and 
inverses  of  these  variables.  It  must  be  admitted  that  some 
pretty  wild  f'unctions  were  generated,  as  will  be  evident  from 
the  actual  fit  developed  which  will  be  presented  below. 

Table  II  lists  the  important  parameters  and  modulii  used  to 
generate  the  113  potential  independent  variables  used. 

The  simplest,  two-term  fit  to  the  strength  data  is; 

S = 0.08278  + 0.64020 

9 

Table  III  shows  the  most  complex,  linear,  least-squares  fit 
which  could  be  constructed.  It  contains  33  terms.  Insofar 
as  the  standard  deviation  is  a good  indication  of  the  quality 
of  a fit,  the  ^fits  obtained  in  this  effort  are  quite  poor.  At 
the  limit  of  statistical  significance,  the  standard  deviation 
could  be  reduced  from  0.29  to  0.20  Kb,  only  about  30  percent 
as  shown  in  Table  IV. 

The  quality  and  characteristics  of  the  fits  are  in- 
dicated in  the  following  figures.  Figure  8 plots  the  measured 
strength  against  the  simple,  two-term  fit  shown  above.  Note 
that  that  fit  does  not  explicitly  involve  the  water  content 
or  saturation.  The  influence  of  saturation  is  shown  in 
Figures  9 and  10  and  water  content  is  examined  in  Figures  11 
and  12.  Finally,  Figure  13  shows  the  best  fit  found  to 
this  strength  data.  It  does  not  differ  noticeably  from  the 
simple  fit  of  Figure  8. 
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Table  IV 

Standard  Deviation  of  Strength  Data 
From  Several  Specific  Fits  to  it 


Simple  average  S = 0.45251  Kb 

a = 0.2921 

Simple  fit  (2  terms) 

a = 0.2159 

Best  fit  (33  terms) 

a = 0.1970 
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OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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SIMPLE  FIT  OF  ALL  DATA 


Figure  8 


OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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Figure  9. 
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OBSERVED  ROCK  STRENGTH 


OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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SIMPLE  FIT  OF  DATA  WITH  MASS  FRACTION  OF  H20  GREATER  THAN  .23 


Figure  11. 
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OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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Figure  12. 
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OBSERVED  ROCK  STRENGTH 


ROCK  STRENGTH  IN  RELATION  TO  OTHER  RHEOLOGICAL  PROPERTIES 
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Figure  13 


3.  CONCLUSIONS 


This  analysis  of  area  12  material  properties  and 
measured  strengths  indicates  that  the  strength  as  repre- 
sented by  stress  difference  at  4 Kb  in  uniaxial  strain  tests 
is  not  strongly  correlated  with  any  other  measured  material 
property.  In  particular,  the  correlation  with  water  content 
is  weak.  One  can  conclude  from  this  work  that  if  strength 
values  are  desired  for  some  purpose,  they  must  be  measured 
or  deduced  by  an  independent  technique.  They  cannot  be 
estimated  reliably  for  tuff  on  the  basis  of  other,  conven- 
tionally-determined properties. 

A strong  caveat  is  certainly  appropriate.  We  all 
recognize  that  the  statement  that  the  dependence  of  material 
strength  on  water  content  is  weak  must  be  viewed  in  context. 
Without  question,  if  one  had  a tuff  in  which  the  water  con- 
tent exceeded,  say,  50  percent,  it  would  be  a very  strange 
material  indeed  and  would  probably  have  no  strength  at  all. 
Four  hundred  and  seventy-one  independent  sets  of  material 
property  data  on  real  materials  were  used  in  this  study. 

For  those  real  materials  an  effort  was  made  to  find  a good 
relation  between  this  definition  of  strength  and  other  par- 
ameters. None  was  found.  That  much  is  true.  One  should  be 
careful  in  drawing  further  conclusions  from  t.his  work. 

Another  class  of  questions  should  be  addressed  in 
the  context  of  material  properties,  residual  stresses  and 
containment. 

(1)  How  meaningful  is  stress  difference  at  4 kbar 
from  uniaxial  strain  tests  as  a fundamental  measure  of 
material  strength?  The  uniaxial  constraint  is  appropriate 
for  a shock  process,  but  not  for  unloading  phenomena.  The 
4 kbar  value  is  of  interest  near  the  shot  point,  but  re- 
sidual stresses  involve  stress  states  closer  to  1.0  kbar. 
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Hopefully  an  answer  to  this  basic  question  will  be  found  in 
the  future. 

(2)  Assuming  this  stress  difference  is  the  value  of 
interest,  is  it  reliably  measured?  This  question  involves 
aspects  of  intrinsic  in-situ  material  variability  and  experi- 
mental technique.  The  question  is  asked  because  a monu- 
mental effort  on  Mighty  Epic  apparently  led  us  to  a wrong 
qualitative  and  quantitative  understanding  of  the  strength 
of  that  material.  Again,  this  question  deserves  our  attention. 
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